Systems and methods to enhance spatial diversity in distributed-input distributed-output wireless systems

ABSTRACT

Systems and methods are described for enhancing the channel spatial diversity in a multiple antenna system (MAS) with multi-user (MU) transmissions (“MU-MAS”), by exploiting channel selectivity indicators. The proposed methods are: i) antenna selection; ii) user selection; iii) transmit power balancing. All three methods, or any combination of those, are shown to provide significant performance gains in DIDO systems in practical propagation scenarios.

RELATED APPLICATIONS

This application is a continuation of the following co-pending U.S. patent application Ser. No. 13/475,598, filed May 18, 2012, which is a continuation-in-part of the following co-pending U.S. patent applications and issued patents:

U.S. application Ser. No. 13/464,648, entitled “System and Methods to Compensate for Doppler Effects in Distributed-Input Distributed Output Systems.”

U.S. application Ser. No. 12/917,257, entitled “Systems And Methods To Coordinate Transmissions In Distributed Wireless Systems Via User Clustering”

U.S. application Ser. No. 12/802,988, entitled “Interference Management, Handoff, Power Control And Link Adaptation In Distributed-Input Distributed-Output (DIDO) Communication Systems”

U.S. Pat. No. 8,170,081, issued May 1, 2012, entitled “System And Method For Adjusting DIDO Interference Cancellation Based On Signal Strength Measurements”

U.S. application Ser. No. 12/802,974, entitled “System And Method For Managing Inter-Cluster Handoff Of Clients Which Traverse Multiple DIDO Clusters”

U.S. application Ser. No. 12/802,989, entitled “System And Method For Managing Handoff Of A Client Between Different Distributed-Input-Distributed-Output (DIDO) Networks Based On Detected Velocity Of The Client”

U.S. application Ser. No. 12/802,958, entitled “System And Method For Power Control And Antenna Grouping In A Distributed-Input-Distributed-Output (DIDO) Network”

U.S. application Ser. No. 12/802,975, entitled “System And Method For Link adaptation In DIDO Multicarrier Systems”

U.S. application Ser. No. 12/802,938, entitled “System And Method For DIDO Precoding Interpolation In Multicarrier Systems”

U.S. application Ser. No. 12/630,627, entitled “System and Method For Distributed Antenna Wireless Communications”

U.S. Pat. No. 7,599,420, issued Oct. 6, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 7,633,994, issued Dec. 15, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 7,636,381, issued Dec. 22, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 8,160,121, issued Apr. 17, 2012, entitled, “System and Method For Distributed Input-Distributed Output Wireless Communications”;

U.S. application Ser. No. 11/256,478, entitled “System and Method For Spatial-Multiplexed Tropospheric Scatter Communications”;

U.S. Pat. No. 7,418,053, issued Aug. 26, 2008, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. application Ser. No. 10/817,731, entitled “System and Method For Enhancing Near Vertical Incidence Skywave (“NVIS”) Communication Using Space-Time Coding.”

BACKGROUND

Prior art multi-user wireless systems add complexity and introduce limitations to wireless networks which result in a situation where a given user's experience (e.g. available bandwidth, latency, predictability, reliability) is impacted by the utilization of the spectrum by other users in the area. Given the increasing demands for aggregate bandwidth within wireless spectrum shared by multiple users, and the increasing growth of applications that can rely upon multi-user wireless network reliability, predictability and low latency for a given user, it is apparent that prior art multi-user wireless technology suffers from many limitations. Indeed, with the limited availability of spectrum suitable for particular types of wireless communications (e.g. at wavelengths that are efficient in penetrating building walls), prior art wireless techniques will be insufficient to meet the increasing demands for bandwidth that is reliable, predictable and low-latency.

BRIEF DESCRIPTION OF THE DRAWINGS

The patent or application file contains at least one drawing executed in color. Copies of this patent or patent publication with color drawing(s) will be provided by the U.S. Patent and Trademark Office upon request and payment of the necessary fee.

A better understanding of the present invention can be obtained from the following detailed description in conjunction with the drawings, in which:

FIG. 1 illustrates one embodiment of a multi-user (MU) multiple antenna system (MAS), or MU-MAS, consisting of a precoding transformation unit.

FIG. 2 illustrates one embodiment in which base transceiver stations (BTSs) are directly connected to a centralized processor (CP).

FIG. 3 illustrates an embodiment in which base transceiver stations (BTSs) are connected via a base station network.

FIG. 4 illustrates another embodiment in which base transceiver stations (BTSs) are connected via a network.

FIG. 5 illustrates simulated transmit covariance matrices for DIDO 6x6 systems with three models.

FIG. 6 illustrates the cumulative density function (CDF) of the diversity metric as well as the symbol error rate (SER) performance as a function of the signal-to-noise ratio (SNR) for three channel models.

FIG. 7 illustrates an exemplary distribution of base transceiver stations.

FIG. 8 illustrates space selectivity indicator (SSI) as a function of the maximum auto-correlation and cross-correlation coefficients of the spatial covariance matrix.

FIGS. 9A-9C illustrate exemplary SSI results for three channel models.

FIG. 10 illustrates a comparison of the cumulative density function (CDF) of the SSI in the three scenarios above.

FIG. 11 illustrates the SNDR of client devices in a DIDO 2x2 system for one measurement set.

FIG. 12 illustrates combined plots containing SNDR, TSI, and SSI.

FIG. 13 illustrates the results from a second channel scenario where RX2 moves from a location with high power imbalance to another with low power imbalance.

FIG. 14 illustrates that even the performance of RX1 (the stationary client) improves as the SSI increases.

FIGS. 15A-15B illustrate the SNDR versus the average SSI (a) and standard deviation of the SSI (b).

FIG. 16 illustrates how SNDR decreases as the standard deviation of the TSI due to deep-fade in time caused by a client's mobility and Doppler effect.

FIG. 17A illustrates the CN defined in equation (7) plotted as a function of the minimum auto-correlation coefficient and maximum cross-correlation coefficient, where every dot is the result of 100 msec of data.

FIG. 17B illustrates the CDF of the SSI defined as the CN.

FIG. 18 illustrates a three dimensional version of FIG. 17 a.

FIG. 19 illustrates the average SNDR as a function of the average CN.

FIGS. 20A-20B illustrate the performance of different order DIDO systems in i.i.d. channels assuming same per-antenna transmit power.

FIGS. 21A-21B illustrate the simulated performance of DIDO 4x4 in i.i.d. channels.

FIGS. 22A-22B illustrate the gain in SNR as a function of the gain of the diversity metric in (9).

FIGS. 23A-23B illustrate the traces of SNDR, SSI (indicated as λ_(min)) and TSI (being the absolute value of the complex channel gain from every transmitter indicated as TX1, . . . , M).

FIGS. 24A-24B illustrate the results obtained over the same route as in FIG. 23, but using two extra antennas.

FIGS. 25A-25B, 26A-26B, 27A-27B, and 28A-28B show the simulated SER performance and CDF of the SSI for DIDO 4x2, 6x4, 8x6 and 10x8 systems.

FIGS. 29A-29B illustrates the average SNDR (over periods of 5 seconds) versus the average SSI.

FIGS. 30A-30B illustrates results for DIDO 4x4 and 6x4.

FIGS. 31A-31B illustrates the average SNDR as a function of the standard deviation of the TSI for DIDO 2x2 and 4x2.

FIGS. 32A-32B illustrates results for DIDO 4x4 and 6x4.

FIG. 33 compares the CDF of the instantaneous SSI for four DIDO systems.

FIG. 34 illustrates the results from one particular measurement set with approximately 20 seconds of data.

FIG. 35 illustrates one embodiment where the SSI target is raised to −15 dB and as a result the SER performance improves.

FIG. 36 illustrates improvements achieved by raising the target SSI.

FIG. 37 illustrates an exemplary scenario where multiple UEs are clustered around one BTS.

FIG. 38 illustrates SER performance of different order DIDO systems for different values of the maximum auto-correlation coefficient.

FIG. 39 illustrates the CDF of the SSI for different values of maximum auto-correlation.

FIG. 40 illustrates one embodiment of a method for balancing the transmit power across all BTSs in the MU-MAS or DIDO system.

FIG. 41 illustrates another embodiment of a method for balancing the transmit power across all BTSs in the MU-MAS or DIDO system.

FIG. 42 illustrates the performance of the transmit power balancing methods in practical outdoor propagation scenarios.

FIG. 43 illustrates the distribution of the condition number with and without power imbalance.

FIGS. 44-46 illustrate the channel traces (SNDR, SSI and TSI) for three different cases: i) DIDO 2x2 without transmit power balancing method; ii) DIDO 2x2 with transmit power balancing method; iii) DIDO 4x2 with transmit power balancing method in combination with antenna selection method.

FIG. 47 illustrates a scenario where a particular BTS is the source of transmit power imbalance.

FIG. 48 illustrates the condition number distribution with and without imbalance.

FIGS. 49-51 depict channel traces for different algorithms.

FIG. 52 illustrates the SER and CDF of the SSI for DIDO systems with 4 clients.

FIG. 53 shows the statistics of the BTSs that have been selected for transmission and their relative usage.

DETAILED DESCRIPTION

One solution to overcome many of the above prior art limitations is an embodiment of Distributed-Input Distributed-Output (DIDO) technology. DIDO technology is described in the following patents and patent applications, all of which are assigned the assignee of the present patent and are incorporated by reference. These patents and applications are sometimes referred to collectively herein as the “related patents and applications.”

U.S. application Ser. No. 13/464,648, entitled “System and Methods to Compensate for Doppler Effects in Distributed-Input Distributed Output Systems.”

U.S. application Ser. No. 12/917,257, entitled “Systems And Methods To Coordinate Transmissions In Distributed Wireless Systems Via User Clustering”

U.S. application Ser. No. 12/802,988, entitled “Interference Management, Handoff, Power Control And Link Adaptation In Distributed-Input Distributed-Output (DIDO) Communication Systems”

U.S. Pat. No. 8,170,081, issued May 1, 2012, entitled “System And Method For Adjusting DIDO Interference Cancellation Based On Signal Strength Measurements”

U.S. application Ser. No. 12/802,974, entitled “System And Method For Managing Inter-Cluster Handoff Of Clients Which Traverse Multiple DIDO Clusters”

U.S. application Ser. No. 12/802,989, entitled “System And Method For Managing Handoff Of A Client Between Different Distributed-Input-Distributed-Output (DIDO) Networks Based On Detected Velocity Of The Client”

U.S. application Ser. No. 12/802,958, entitled “System And Method For Power Control And Antenna Grouping In A Distributed-Input-Distributed-Output (DIDO) Network”

U.S. application Ser. No. 12/802,975, entitled “System And Method For Link adaptation In DIDO Multicarrier Systems”

U.S. application Ser. No. 12/802,938, entitled “System And Method For DIDO Precoding Interpolation In Multicarrier Systems”

U.S. application Ser. No. 12/630,627, entitled “System and Method For Distributed Antenna Wireless Communications”

U.S. Pat. No. 7,599,420, issued Oct. 6, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 7,633,994, issued Dec. 15, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 7,636,381, issued Dec. 22, 2009, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. Pat. No. 8,160,121, issued Apr. 17, 2012, entitled, “System and Method For Distributed Input-Distributed Output Wireless Communications”;

U.S. application Ser. No. 11/256,478, entitled “System and Method For Spatial-Multiplexed Tropospheric Scatter Communications”;

U.S. Pat. No. 7,418,053, issued Aug. 26, 2008, entitled “System and Method for Distributed Input Distributed Output Wireless Communication”;

U.S. application Ser. No. 10/817,731, entitled “System and Method For Enhancing Near Vertical Incidence Skywave (“NVIS”) Communication Using Space-Time Coding.”

To reduce the size and complexity of the present patent application, the disclosure of some of the related patents and applications is not explicitly set forth below. Please see the related patents and applications for a full detailed description of the disclosure.

1. SYSTEM MODEL

Described below is a multi-user (MU) multiple antenna system (MAS), or MU-MAS, consisting of a precoding transformation unit 101, a network 102 and M transceiver stations 103 communicating wirelessly to N client devices UE1-UE4, as depicted in FIG. 1. The precoding transformation unit 101 receives N streams of information with different network contents (e.g., videos, web-pages, video games, text, voice, etc., streamed from Web servers or other network sources C1-C5) intended for different client devices. Hereafter, we use the term “stream of information” to refer to any stream of data sent over the network containing information that can be demodulated or decoded as a standalone stream, according to certain modulation/coding scheme or protocol, to produce certain voice, data or video content. In one embodiment, the stream of information is a sequence of bits carrying network content that can be demodulated or decoded as a standalone stream. In one embodiment, this network content is delivered to the precoding transformation unit 101 via a network. Any type of network access technology may be used including wireline and wireless. Additionally, the network may be a local area network (e.g., LAN, WLAN, etc.), wide area network, the Internet, or any combination thereof.

In one embodiment, the precoding transformation unit 101 processes the channel state information (CSI) for each communication channel established with each client device UE1-UE4 to produce a precoding transformation. In another embodiment, channel quality information (e.g., signal-to-noise ratio, etc) or statistical channel information (e.g., spatial covariance matrix, etc.) are used to compute the precoding transformation. The precoding transformation can be linear (e.g., zero-forcing [1], block-diagonalization [2], matrix inversion, etc.) or non-linear (e.g., dirty-paper coding [3-5] or Tomlinson-Harashima precoding [6-7]).

In one embodiment, the precoding transformation unit 101 utilizes the precoding transformation to combine (according to certain algorithm) the N streams of information from the network content C1-C5 into M streams of bits. Hereafter, we use the term “stream of bits” to refer to any sequence of bits that does not necessarily contain any useful bit of information and as such cannot be demodulated or decoded as a standalone stream to retrieve the network content. In one embodiment of the invention, the stream of bits is the complex baseband signal produced by the precoding transformation unit and quantized over given number of bits to be sent to one of the M transceiver stations 103. In one embodiment, the M streams of bits are sent from the precoding transformation unit to the M transceiver stations 103 via the network 102 (which may be a wireline/wireless, Internet, wide area network, or local area network, or any combination thereof).

Finally, the M transceiver stations 103 send the streams of bits to the client devices UE1-UE4 that recover the streams of information and demodulate the network content. Note that the number of clients K in the system can be any value. For example, if K>M the extra (K−M) clients are multiplexed via different techniques described in the related patents and applications and in the prior art (e.g., TDMA, FDMA, OFDM, CDMA, etc.). Also, if K<=M but K<N, more than one stream of information is available for some of the client devices. Those client devices can demodulate multiple streams of information if they are equipped with multiple antennas by using existing MIMO or DIDO techniques.

One important feature of the present invention is that the MU-MAS transforms the streams of information into streams of bits sent over the network to the transceiver stations 103, such that the client devices UE1-UE4 can recover the stream of information when receiving the streams of bits simultaneously from all transceiver stations. We observe that, unlike prior art, the M streams of bits sent through the network are combinations of some or all N streams of information. As such, if a client device had to receive the stream of bits from only one of the M transceiver stations (even assuming good link quality and SNR from that station to the client), that information would be completely useless and it would be impossible to recover the original network content. It is only by receiving the streams of bits from all or a subset of the M transceiver stations that every client device can recover the streams of information and demodulate the network contents C1-C5.

In one embodiment of the invention, the MU-MAS is a distributed-input distributed-output (DIDO) system consisting of a centralized processor (CP) 201, base transceiver stations (BTSs) 203, and user equipment (UEs) UE1-UE4 as shown in FIG. 2 and described in the related patents and applications referenced above. The BTSs can either be directly connected to the CP 201 as in FIG. 2 or via the base station network (BSN) 301 as depicted in FIG. 3. In another embodiment, the network contents C1-C5 and the BTSs 203 are both connected to the CP 201 via the same network 401 as in FIG. 4, which may be a wireline/wireless local area network, wide area network, and/or the Internet.

For client devices to reliably recover the network content from the received streams of information, the wireless channel must have a sufficient number of degrees of freedom or equivalently must have high spatial diversity. Spatial diversity depends on the distribution in space of the transceiver stations 203 and the client devices UE1-UE4 as well as the spatial distribution of multi-paths in the propagation environment (or channel angular spread). Described below are different metrics to evaluate the spatial diversity of the wireless channel that will be used in the techniques and methods described later on in the present application.

2. DIVERSITY METRICS AND CHANNEL MODELS

The received signal at target client k is given by

$\begin{matrix} {r_{k} = {{H_{k}W_{k}s_{k}} + {H_{k}{\sum\limits_{\underset{u\; \neq \; k}{u = 1}}^{K}{W_{u}s_{u}}}} + n_{k}}} & (1) \end{matrix}$

where k=1, . . . , K, with K being the number of clients. Moreover, r_(k)εC^(R×M) is the vector containing the receive data streams at client k, assuming M transmit DIDO antennas and R receive antennas at the client devices; s_(k)εC^(N×1) is the vector of transmit data streams to client k in the main DIDO cluster; s_(u)εC^(N×1) is the vector of transmit data streams to client u in the main DIDO cluster; n_(k)εC^(N×1) is the vector of additive white Gaussian noise (AWGN) at the R receive antennas of client k; H_(k)εC^(R×M) is the DIDO channel matrix from the M transmit DIDO antennas to the R receive antennas at client k; W_(k)εC^(M×R) is the matrix of DIDO precoding weights to client k in the main DIDO cluster; W_(u)εC^(M×R) is the matrix of DIDO precoding weights to client u in the main DIDO cluster.

To simplify the notation without loss of generality, we assume all clients are equipped with R receive antennas and that there are M DIDO distributed antennas with M≧(R·K). If M is larger than the total number of receive antennas, the extra transmit antennas are used to pre-cancel interference to the target clients or to improve link robustness to the clients within the same cluster via diversity schemes described in the related patents and applications, including U.S. Pat. Nos. 7,599,420; 7,633,994; 7,636,381; and application Ser. No. 12/143,503.

The DIDO precoding weights are computed to pre-cancel inter-client interference. For example, block diagonalization (BD) precoding described in the related patents and applications, including U.S. Pat. Nos. 7,599,420; 7,633,994; 7,636,381; and application Ser. No. 12/143,503 and [2] can be used to remove inter-client interference, such that the following condition is satisfied in the main cluster

H _(k) W _(u)=0^(R×R) ;∀u=1, . . . ,K;withu≠k.  (2)

Substituting conditions (2) into (1), we obtain the received data streams for target client k, where inter-user interference is removed

r _(k) =H _(k) W _(k) s _(k) +n _(k).  (3)

We define the effective channel matrix of user k as

{tilde over (H)} _(k) =H _(k) W _(k).  (4)

One embodiment of the invention defines the diversity metric as the minimum over all clients of the minimum singular values of the effective channel matrices in (4)

$\begin{matrix} {\lambda_{\min} = {\min\limits_{{k\; = \; 1},\; \ldots \;,\; K}{{\lambda_{\min}^{(k)}\left( {\overset{\sim}{H}}_{k} \right)}.}}} & (5) \end{matrix}$

Another embodiment uses the minimum or maximum singular value or the condition number of the composite DIDO channel matrix obtained by staking the channel matrices from every client as

$\begin{matrix} {\overset{\_}{H} = {\begin{bmatrix} H_{1} \\ \vdots \\ H_{K} \end{bmatrix}.}} & (6) \end{matrix}$

The condition number (CN) is defined as the ratio between the maximum and the minimum singular value of the composite DIDO channel matrix as

$\begin{matrix} {{C\; N} = {\frac{\lambda_{\max}\left( \overset{\_}{H} \right)}{\lambda_{\min}\left( \overset{\_}{H} \right)}.}} & (7) \end{matrix}$

Next, we define different channel models that will be used to simulate the performance of the system and methods described in this application in realistic propagation conditions. We employ the well known Kronecker structure [8,9] and model the spatial covariance matrix with cross-correlation entries r_(i,j)=ρ_(c) ^(|i−j|) with i≠j, and auto-correlation entries given by

$\begin{matrix} {r_{i,i} = \left\{ {\begin{matrix} \rho_{a} & {i = 1} \\ \frac{\left( {N - \rho_{a}} \right)\left( {N - i + 1} \right)}{\sum\limits_{n = 1}^{N - 1}n} & {i > 1} \end{matrix}.} \right.} & (8) \end{matrix}$

The advantage of the model in (8) is that it allows us to write the diagonal entries of the spatial covariance matrix as a function of only one parameter ρ_(a).

We define three different channel models for the transmit spatial covariance matrix: i) “i.i.d. model” with ρ_(c)=0.00001, ρ_(a)=1 that approximates the independent identically distributed models; ii) “high cross-correlation model” with ρ_(c)=0.8, ρ_(a)=1 to simulate wireless systems where the antennas have equal transmit power and are in close proximity to each other (e.g., corner case in MIMO systems) thereby yielding high cross-correlation coefficients; iii) “high auto-correlation model” with ρ_(c)=0.00001, ρ_(a)=5.9 to simulate wireless systems with antennas distributed over a large area to yield low spatial correlation, but with one antenna overpowering all the others due to its close proximity to all clients (e.g., corner case in DIDO systems). Simulated transmit covariance matrices for DIDO 6x6 systems with these three models are shown in FIG. 5. In all the results presented hereafter, we assume the receive covariance matrix is identity, since the clients are assumed to be spread over large area, several wavelengths apart from one another.

FIG. 6 shows the cumulative density function (CDF) of the diversity metric (i.e., minimum singular value) as well as the symbol error rate (SER) performance as a function of the signal-to-noise ratio (SNR) for the three channel models described above. We observe the SER performance in the “high cross-correlation model” and “high auto-correlation model” degrades due to lack of spatial diversity. In the “high cross-correlation model”, lack of diversity is due to high spatial correlation across the transmit antennas of the MIMO array. In the “high auto-correlation model”, reduced diversity is due to transmit power imbalance on one of the transmit antennas over the others. One way to improve spatial diversity in the “high cross-correlation model” is to space antennas far apart, which can be prohibitive in practical MIMO systems. In the “high auto-correlation model” diversity can be increased by balancing transmit power across the distributed antennas as described later in this patent application.

We collected the diversity metric in DIDO systems for a variety of propagation conditions. In the experimental campaign, we used the DIDO BTSs installed in different buildings in downtown Palo Alto, as shown in FIG. 7. We began by measuring the “space selectivity indicator” (SSI) in (5) for DIDO 2x2 systems in a variety of propagation scenarios and with different combinations of transmit BTSs and receive antenna locations around downtown Palo Alto.

FIG. 8 shows the SSI as a function of the maximum auto-correlation and cross-correlation coefficients of the spatial covariance matrix. Each dot is obtained by averaging the spatial covariance matrix over a period of 5 seconds, which is enough to average out the fading effect at the speed considered in the experiments (i.e., stationary clients as well as clients moving at 3 mph). We observe the highest values of SSI (that indicate high channel spatial selectivity) are obtained when the maximum cross-correlation is “0” (i.e., low spatial correlation between transmit antennas, due to large physical spacing between antennas or high channel angular spread) and the maximum auto-correlation is “1” (i.e., good power balance across transmit antennas). Any scenario that deviates from these two cases yields low values of SSI and low spatial diversity.

FIGS. 9a-9b shows the SSI results for the three channel models described above. The “i.i.d.” case in FIG. 9a indicates scenarios where the DIDO BTSs were physically faraway from one another; the “high cross-correlation” case was obtained by spacing the transmit antennas one wavelength apart while allowing the clients to move around anywhere in the coverage area; the “high auto-correlation” case was obtained by placing all clients in proximity of antenna 10 in FIG. 7, such that it would overpower the other transmitter. The pink circle in every plot indicates the average of the auto- and cross-correlation coefficients. In FIG. 9b the average cross-correlation increases from FIG. 9a due to small antenna spacing and the average auto-correlation decreases due to transmit power balance. Vice versa, in FIG. 9c the average auto-correlation increases due to transmit power imbalance and the average cross-correlation decreases due to larger antennas spacing between BTSs.

FIG. 10 compares the cumulative density function (CDF) of the SSI in the three scenarios above. We observe that the i.i.d. case yields the best performance over the other two channel conditions.

Next, we analyze how the signal-to-noise-plus-distortion ratio (SNDR) and SER performance of DIDO systems varies as a function of spatial and temporal variations. The spatial variations are measured via the above defined SSI. The temporal variations are measured through the “time selectivity indicator” (TSI). One embodiment of the inventions defines the TSI as the absolute value of the sum of the complex channel gain from some or all transmit antennas in the DIDO system. Any other metric tracking channel variations, deep-fade rate or duration can be used as TSI. The top row in FIG. 11 shows the SNDR of both client devices in the DIDO 2x2 system for one measurement set. The second row depicts the TSI: we observe the TSI of client 1 is flat because it is stationary, whereas the TSI of client 2 fades over time due to client's mobility. The third row shows the SSI as well as the auto-correlation coefficients for each of the two transmitters. We observe that when the two auto-correlation coefficients are close to “1”, the SSI increases.

FIG. 12 combines the above results all on one plot. For the first receiver (i.e., RX1) we observe the SNDR trace fades even though the client is stationary and the TSI is flat. In fact, due to the mobility of RX2, the SSI varies over time and those variations produce SNDR fades also for RX1. Moreover, we observe that the SSI may fade independently on the TSI. In fact TSI fades indicate poor signal quality from both transmitters due to destructive interference of multipaths, but the DIDO channel may still have enough spatial degrees of freedom (i.e., large SSI) to support multiple spatial data streams. Alternatively, SSI fades indicate the DIDO channel matrix is singular and cannot support multiple parallel data streams, but the signal quality from all transmit antennas may still be good, yielding large TSI. The present invention uses TSI, SSI, or a combination of the two metrics to evaluate the channel diversity and adjust system parameters to enhance diversity.

The SSI can be used to measure and predict the areas of coherence in DIDO systems. For example, one embodiment of the invention measures the SSI, keeps track of it over time, and predicts its future behavior. Based on that prediction, it adapts both transmit and receive system parameters (e.g., number of BTSs to employ for transmission or number of client devices to receive data streams).

FIG. 13 shows the results from a second channel scenario where RX2 moves from a location with high power imbalance to another with low power imbalance. The result of this variation is SSI improvement. FIG. 14 shows the summary for the same measurement scenario: we observe even the performance of RX1 (the stationary client) improves as the SSI increases.

We compared the SNDR performance against the SSI in a large set of propagation conditions. FIGS. 15a-15b depicts the SNDR versus the average SSI in a) and standard deviation of the SSI in b). Each dot represents one measurement collected over a period of 5 seconds. The red solid line represents the linear regression of all measurement points, where the upper and lower 10% of the measurements has been removed from the regression calculation. We observe the average SNDR increases as a function of the average SSI in FIG. 15a due to large spatial diversity available in the wireless channel. For example, scenarios characterized by large power imbalance yield low average SSI resulting in low SNDR. Moreover, in FIG. 15b the average SNDR decreases as a function of the standard deviation of the SSI due to frequent deep-fades of the SSI due to client mobility. Note that in practical systems the average and standard deviation of the SSI can be computed via running average techniques or methods using forgetting factor for efficient memory use and reduction in computational complexity.

FIG. 16 shows similar results for the temporal channel variations. The SNDR decreases as the standard deviation of the TSI due to deep-fade in time caused by the client's mobility and Doppler effect.

Another embodiment of the invention uses the condition number (CN) as SSI. The CN defined in equation (7) is plotted as a function of the minimum auto-correlation coefficient and maximum cross-correlation coefficient in FIG. 17b , where every dot is the result of 100 msec of data. Contrarily to the minimum singular value of the effective channel matrix, lower CN indicates channels with high spatial diversity. FIG. 17a shows the CDF of the SSI defined as the CN. FIG. 18 depicts the three dimensional version of FIG. 17 a.

FIG. 19 shows the average SNDR as a function of the average CN. Every dot represents an average over 20 seconds of data. We observe that the SNDR degrades as the value of the average CN increases.

3. METHODS FOR TRANSMIT ANTENNA SELECTION

One way to increase the spatial degrees of freedom in a wireless link is to add more transmit antennas than the number of clients in the system and select among the antennas that satisfy a certain SSI performance target. This algorithm is known as transmit antenna selection as described in [10] and our previous patent application U.S. Pat. No. 7,636,381. In one embodiment, all possible combination of transmit antenna subsets are first identified. Then the SSI is computed for each of the antenna sets. Finally the set that maximizes the diversity metric or SSI is chosen as optimal transmit antenna subset.

FIGS. 20a-20b show the performance of different order DIDO systems in i.i.d. channels assuming same per-antenna transmit power. The SSI degrades for increasing number of transmit antennas as the CDF shifts to the left going from 2x2 to 8x8, but the SER performance is similar for any order DIDO.

FIGS. 21a-21b show the simulated performance of DIDO 4x4 in i.i.d. channels. Antenna selection provides significant gains in SNR depending on the target SER. For example, at SER target of 1% by adding two extra antennas the gain is 12 dB or at a target of 0.1% that gain increases up to 18 dB. Also, FIG. 21b shows that the CDF of the SSI improves with antenna selection due to enhanced spatial diversity. Note that in FIG. 21 we plot the maximum (over all possible transmit antenna subsets) of the minimum singular value of the effective channel matrix. We define the mean value of the CDF in FIG. 21b in decibels as

D _(dB)=20 log₁₀(E{max(λ_(min))})  (9)

FIGS. 22a-22b show the gain in SNR as a function of the gain of the diversity metric in (9). The values in the table a) are obtained from the simulated SER performance in FIG. 21a . In FIG. 22b we observe close to linear relation between the two gains. In one embodiment of the invention, the average SSI is used to decide whether to employ the selection algorithm or not. In fact, antenna selection algorithms require additional computational complexity as the SSI must be computed over all antenna subsets. Understanding under what channel conditions the antenna selection algorithm is really needed, allows to turn off the algorithm when unnecessary, thereby improving computational complexity of the system. For example, if the average SSI is above a certain threshold, there is no need to trigger the antenna selection algorithm and a fixed number of antennas are used for transmission. If the average SSI decreases (e.g., due to client's mobility) the algorithm will trigger the antenna selection algorithm in such a way that the best antenna subset is chosen for transmission over the wireless link.

The SSI threshold can be pre-calculated by analyzing experimental data from practical measurements. For example, FIGS. 23a-23b show the traces of SNDR, SSI (indicated as A_(min)) and TSI (being the absolute value of the complex channel gain from every transmitter indicated as TX1, . . . , M). Note that we intentionally subtracted 20 dB from the TSI traces to fit all traces into the same plot while avoiding overlaps. In this experiment, the first client RX1 is stationary, whereas the second RX2 is mobile. We observe that even for the stationary client the SNDR trace varies over time due to fades in the SSI trace. In particular, every time the SSI falls below −10 dB, the SNDR undergoes deep-fades. We choose −10 dB as a threshold for the SSI. This invention is not limited to this value of SSI and other values may be chosen based on different performance criteria. For the mobile client, deep-fades are caused by either SSI fades or TSI fades. As observed before, these two types of fades are uncorrelated and may occur at different times.

FIGS. 24a-24b show the results obtained over the same route as in FIGS. 23a-23b , but using two extra antennas. Note that the SNDR trace is not aligned with the one in FIGS. 24a-24b because the instantaneous channel varies from one experiment to the next due to fast-fading effects. We observe that, by adding two extra antennas and running the transmit antenna selection algorithm, it is possible to remove deep-fades from the SSI trace and improve SNDR performance of both clients. FIG. 24a shows that the stationary client does not undergo any SNDR deep fade. FIG. 24b shows that the SNDR fades of the mobile client are only due to TSI, whereas the SSI fades are completely removed.

One embodiment of the invention scans through the available transmit antenna subsets until the first one that provides SSI above the predefined threshold is reached. Once that subset is found, the search stops thereby reducing the computational complexity of the algorithm.

In FIGS. 23a-23b and 24a-24b we observed that the SSI exhibits a structured behavior with periods of deep-fade that alternate to periods of high gain. This temporal correlation between consecutive samples of SSI can be exploited to reduce the complexity of the antenna subset selection algorithm. In one embodiment, the same antenna subset is preserved until the SSI undergoes a deep-fade and its value drops below the predefined threshold. In another embodiment, the system keeps track of the statistics of the antenna subset that have been chosen in the past and only selects those for future transmissions.

Another way to reduce computational complexity of the system is to reduce the number of combinations of transmit antennas to be chosen across with the antenna selection method. FIGS. 25a-25b, 26a-26b, 27a-27b and 28a-28b show the simulated SER performance and CDF of the SSI for DIDO 4x2, 6x4, 8x6 and 10x8 systems. All these systems employ two extra antennas than the number of clients. Performance is shown for different numbers Ns of antenna subsets. We observe that, for any DIDO order, 10 or less subsets of transmit antennas suffice to approximate closely the SER performance of the same system using all possible combinations of antenna subsets. Reducing the number of antenna subsets can yield a significant reduction in computational complexity as the SSI does not need to be computed over all antenna subsets. One embodiment of the invention selects a limited number of subsets as a means to reduce the computational complexity of the system, while maintaining system performance close to ideal performance.

One embodiment of the invention uses combination of SSI and TSI to select the optimal antenna subset. For example, the antenna subset that provides the maximum SSI and TSI is selected. Another embodiment defines a first selection phase that identifies all antenna subsets that provide SSI above the predefined threshold. Then, a second selection phase chooses the subset that yields the largest TSI. Alternatively, another threshold is defined for the TSI and the subset that satisfies both SSI and TSI thresholds is selected.

All the methods and results described above for single-carrier systems can be directly extended to multi-carrier and/or OFDM systems by defining “frequency selectivity indicator” (FSI). For example, in OFDM systems every tone experiences a frequency flat channel. Then all methods described above can be applied on a tone-by-tone basis. In another embodiment, different combinations of SSI, TSI and FSI are employed to select the optimal antenna subset according to the criteria defined above.

Finally, we show the performance of antenna selection algorithms in a variety of propagation conditions. FIGS. 29a-29b depict the average SNDR (over periods of 5 seconds) versus the average SSI. Large average SSI indicates channels with high spatial diversity, thereby yielding large average SNDR. We observe that two extra antennas used for antenna selection going from DIDO 2x2 to 4x2 increases both average SNDR and SSI. Similar results are shown in FIGS. 30a-30b for DIDO 4x4 and 6x4.

FIGS. 31a-31b depicts the average SNDR as a function of the standard deviation of the TSI for DIDO 2x2 and 4x2. High standard deviation indicates large time selectivity due to Doppler effects that degrade the SNDR performance. Similar results are showed in FIGS. 32a-32b for DIDO 4x4 and 6x4.

Finally, FIG. 33 compares the CDF of the instantaneous SSI for all four DIDO systems considered above. We observe that 4x4 has worse CDF performance than 2x2 due to reduced degrees of freedom when switching to higher order DIDO. In both cases, adding 2 extra antennas with transmit selection algorithms yield significant improvement in SSI performance.

4. METHODS FOR USER SELECTION

In one embodiment, spatial diversity is enhanced in DIDO channels via user selection. In this embodiment, if there are not enough degrees of freedom in the wireless channel for the given number of transmit antennas available in the system, then the system drops transmission to one or multiple clients. This technique may employ the SSI to measure the spatial diversity in the wireless link. When the SSI falls below a predefined threshold, one or multiple clients are dropped.

In one embodiment of the invention, the fastest moving client is dropped. In fact, the client experiencing the highest Doppler effect is most likely to undergo deep-fades. Another embodiment utilizes the TSI and FSI to select the client with lower channel quality and drops that client. When the client is dropped, the bits transmitted over that period are corrupted and those bits can be recovered via forward error correction (FEC) coding. Another embodiment utilizes alternative multiplexing technique such as TDMA, FDMA, OFDMA or CDMA to serve the dropped clients.

FIG. 34 shows the results from one particular measurement set with approximately 20 seconds of data. The first row depicts the measured SNDR trace for the two clients denoted as RX1 (stationary client) and RX2 (mobile client); the second row is the simulated SNDR with the target fixed to 10 dB to demodulate 4-QAM constellations reliably; the third row is the simulated SER; finally the fourth row depicts the SSI and the auto-correlation coefficients. We observe that even if RX1 is stationary, its SNDR drops below the target due to lack of spatial diversity, as indicated by low SSI. If the SSI target is too low (i.e., −20 dB in the figure) the user selection algorithm is never triggered and RX1 undergoes severe SER degradation.

FIG. 35 shows the case where the SSI target is raised to −15 dB and as a result the SER performance improves. Further improvement is achieved by raising the target SSI further up −10 dB as in FIG. 36, in which case the SER for RX1 is reduced to zero throughout the duration of the measurement. In this case the SSI threshold is determined based on the SER performance, but this invention is not limited to that and any other performance criterion can be used for that.

5. METHODS FOR TRANSMIT POWER BALANCING

Transmit power imbalance occurs when most or all of the clients are around one BTS and far from all the others, such that one BTS overpowers the others. Transmit power imbalance reduces channel spatial diversity (i.e., decreases the SSI), thereby adversely affecting system performance. One exemplary scenario is shown in FIG. 37 where multiple UEs 3701 (identified as squares) are clustered around one particular BTS 3702 (identified with a circle) and located far away from the other BTSs. This scenario would happen, for example, when there is an event in one location in which the group of clients are participating, and all other BTSs are far away. One embodiment of the invention adaptively adjusts the power of the BTSs in such a way that the power received at all clients from all BTSs is balanced. In one embodiment of the invention, the power of the BTS that is overpowering all the others is reduced until the power received by the clients balances the power received from all other BTSs. In another embodiment of the invention the power from all other BTSs is increased until the received power level from all BTSs to every client is balanced.

In TDD systems in which channel reciprocity is exploited, the channel state information (CSI) for the downlink is obtained from the uplink. The uplink training signal is quantized by the ADC at the receiver of the BTS and, as such, it has limited dynamic range, depending on the number of bits of the ADC. If all clients are clustered around one of the BTSs, the CSI for that BTS will have a much larger amplitude than the one from all the others and, as such, it will make the DIDO channel matrix singular and limit the spatial degrees of freedom of the link. That is the effect of transmit power imbalance. In FDD systems or TDD systems that do not exploit channel reciprocity, the same issue manifests at the receiver of the client devices also equipped with ADC. Moreover, the CSI may need to be quantized or mapped into bits via limited feedback techniques, before being sent over the wireless link. That quantization again limits the dynamic range for the CSI and yields a power imbalance when one of the BTSs overpowers the other. Embodiment of the invention described herein employ techniques for preventing power imbalance in MU-MAS and DIDO systems.

As shown in FIG. 13, one way to identify transmit power imbalance is by looking at the auto-correlation coefficient ρ_(a): when the auto-correlation value approximates the number of the BTSs (assuming the transmit spatial covariance matrix is normalized with trace equal to the number of BTSs) the system undergoes transmit power imbalance. For example, in a power imbalanced DIDO 4x4 system, one auto-correlation would be close to “4” and all other auto-correlation coefficients would be close to zero. Contrarily, in a perfectly balanced system, all auto-correlation coefficients would be “1”.

Transmit power imbalance adversely affects the performance of the system. For example, FIG. 38 shows the SER performance of different order DIDO systems for different values of the maximum auto-correlation coefficient. As the maximum auto-correlation decreases to “1”, the SER performance approaches the ideal i.i.d. case. These SER results can be used to define thresholds that distinguish balanced systems from imbalanced systems. These auto-correlation thresholds can be determined through numerical, analytical, or empirical methods. For example, in FIG. 38 the thresholds are chosen such that the SER performance does not degrade more than 3 dB from the ideal i.i.d. performance. The invention, however, is not limited to this performance criterion and any other criteria that measure the system performance can be used. Another embodiment of the invention employs a hysteresis loop where two different thresholds are defined for the auto-correlation coefficient as in the table in FIG. 38.

FIG. 39 shows the CDF of the SSI for different values of maximum auto-correlation. We observe that increasing the maximum auto-correlation yields worse SSI performance due to reduced spatial diversity.

Embodiments of the invention propose different methods for balancing the transmit power across all BTSs in the MU-MAS or DIDO system. These methods can be executed at a regular rate. In one embodiment, the proposed methods run every execution cycle. However, depending on the constraints of the system being used, a lower rate may be used. Hereafter, we described these methods in details.

One embodiment of the invention aims to keep the transmit power of each BTS at the maximum possible level, while staying within the auto-correlation thresholds. We define two different thresholds, as shown in FIG. 38. The upper threshold, MAX_AUTO_CORR, represents the point at which the auto correlation number results in a significant drop in spatial diversity. If the auto-correlation number goes above this threshold, there will be a large drop in system performance.

The lower threshold, MIN_AUTO_CORR acts as a buffer to prevent the system from changing power settings too often. If a given BTS has an auto correlation number below MIN_AUTO_CORR, it can safely increase its transmit gain value (assuming transmit gain is not already set to its maximum). Note that the transmit gain may be the analog gain of the power amplifier in the RF chain and/or the digital gain corresponding to a certain level of the DAC. If the auto-correlation is between the MIN_AUTO_CORR and MAX_AUTO_CORR, no action is taken. If the power was to be increased in this instance, it could increase the auto-correlation number until it was above the MAX_AUTO_CORR, at which point the power would be decreased until it was below the MAX_AUTO_CORR, etc. This effect would cause the power to be changing constantly, which is inefficient and may potentially cause performance degradation.

One embodiment of a method is illustrated in FIG. 40 and its associated pseudo-code is described as follows:

BEGIN INITIALIZE txGain for each BTS SET highestAutoCorrNum = 0 SET K = 0 REPEAT WHILE K < number of BTSs IF auto correlation number for BTS K > highestAutoCorrNum THEN SET maxAutoCorrNum = auto correlation number of BTS K SET N = K END IF INCREMENT K END REPEAT IF highestAutoCorrNum > MAX_AUTO_CORR AND txGain for BTS N > MIN_TX_GAIN THEN Decrease the txGain for BTS N by TX_GAIN_STEP Set radio TX gain for BTS N = txGain for BTS N SET K = 0 REPEAT WHILE K < number of BTSs IF txGain for BTS K < MAX_TX_GAIN AND auto correlation number for BTS K < MIN_AUTO_CORR THEN INCREASE txGain for BTS K by TX_GAIN_STEP SET radio TX gain for BTS K to txGain for BTS K END IF INCREMENT K END REPEAT END

In summary, this method first determines which BTS has the highest correlation. That correlation value is saved, along with the index of the corresponding BTS. Then, if the highest correlation is above the upper threshold, the transmit gain is decreased. The transmit gain will not decrease below a defined minimum. Then, for each BTS, the transmit gain is increased if the highest correlation is below the lowest value. If the highest auto-correlation number is between the two thresholds, no action is taken. This is the target mode of operation of the proposed method.

Turning to the specific details of FIG. 40, at 4001, the variables highestAutoCorrNumber and K are initialized to zero. Steps 4002 and 4004 ensure that the loop repeats for each BTS. At 4002, if the current value of K is currently less than the number of BTSs, then at 4003, a determination is made as to whether the autocorrelation number for BTS K is greater than the current highest auto-correlation number. If so, then at 4005 the variable highestAutoCorrNum is set to the auto-correlation number for BTS K (i.e., BTS K has the highest auto-correlation number) and control variable N is set equal to K.

At 4006, if the highestAutoCorrNum is greater than the maximum auto-correlation (MAX_AUTO_CORR) and the transmit gain (txGain) for BTS N is greater than the minimum transmit gain (MIN_TX_GAIN) then, at 4008, the transmit gain for BTS N is decreased using a specified step size (TX_GAIN_STEP) and the txGain of BTS N's radio is set to the new txGain value.

At 4009, the control value K is set equal to zero. Step 4010 ensures that each BTS is addressed by the loop of steps 4011-4012. That is, if K is currently less than the number of BTSs (i.e., if all BTSs have not been analyzed) then, at 4011, a determination is made as to whether the auto-correlation number for BTS K is less than the minimum auto-correlation (MIN_AUTO_CORR) and the txGain for BTS K is less than the maximum allowable transmit gain value (MAX_TX_GAIN). If both conditions are met then, at 4012, the transmit gain for BTS K is increased by the predefined step size (TX_GAIN_STEP) and the new txGain is set on BTS K's radio. The control value K is incremented at 4013 and, at 4010, if K is equal to the number of BTSs (i.e., each BTS has been analyzed), the process terminates.

In another embodiment of the invention, auto-correlation values are mapped to transmit gain values. One embodiment uses a linear mapping, shown below. Although a linear mapping is simple to implement, the adverse effect of the auto-correlation on system performance does not scale linearly. Typically, system performance is significantly affected only after the auto-correlation number reaches some fraction of its maximum value. For example, DIDO 2x2 performance is seriously affected only when the maximum auto-correlation is above 1.95 (or 97.5% of its maximum value). Another mapping algorithm may utilize an exponential function or another power function designed to operate in these ranges, rather than a linear function.

One embodiment of the method is illustrated in FIG. 41 and its pseudo-code is described as follows:

BEGIN INITIALIZE txGain for each BTS SET K = 0 REPEAT WHILE K < number of BTSs SET autoCorr = auto correlation number for BTS K SET txGain for BTS K = (MAX_TX_GAIN − MIN_TX_GAIN) * (1 − autoCorr/nTX) + MIN_TX_GAIN INCREMENT K END REPEAT

This method takes an auto-correlation number and scales it directly into a transmit gain value. Most of the complexity in the method is to allow different orders of DIDO and different values of MIN_TX_GAIN and MAX_TX_GAIN. For example, the simplest form of the equation for a DIDO 2x2 system with transmit gain that ranges between A and B would be:

$\begin{matrix} {\left( {B - A} \right)*\left( {1 - \frac{\rho_{a}}{2}} \right)} & (9) \end{matrix}$

For example, an auto-correlation value of 2 (highest value for DIDO 2x2) would result in the transmit gain for that BTS being set to A=0 (lowest transmit power), while an auto correlation value of 0 (lowest value for DIDO 2x2) would result in the transmit gain for that BTS being set to B=30 (highest transmit power). It should be noted that both of these cases indicated extreme power imbalance. In the first case (ρ_(a)=2.0), this BTS is being received too strongly across the UEs. In the second case (ρ_(a)=0.0), the other BTS is being received too strongly. A perfectly balanced system, with ρ_(a)=1.0 for both BTSs, would result in the transmit gain staying at 15 (being the default value), as desired.

Turning to the specifics of FIG. 41, at 4101, the control variable K is initialized to 0. At 4102, if K is less than the number of BTSs being observed, then the operations set forth in 4103 are performed. Specifically, the variable autoCorr is set equal to the current auto correlation number for BTS K. In addition, the variable txGain for BTS K is set equal to the difference between the maximum transmit gain value and the minimum transmit gain value (MAX_TX_GAIN −MIN_TX_GAIN) multiplied by (1−autoCorr/nTX) and added to the minimum transmit gain value (MIN_TX_GAIN). Control variable K is then incremented until K=the number of BTSs. The process then terminates.

Both of the previous methods are designed to adjust the transmit gain of every BTS within a single step. Another embodiment of the invention defines a method that always adjusts the power of only two BTSs. With this method, however, in certain scenarios one or more of the BTSs could remain at low transmit power setting for long periods of time. Thus, in practical systems this method would be combined with an algorithm similar to Method 1 (using thresholds as in FIG. 40) where the power of each BTS is increased if the auto-correlation number for that BTS is below MIN_AUTO_CORR.

The pseudo-code for Method 3 described above is as follows:

BEGIN INITIALIZE txGain for each BTS SET highestAutoCorrNum = 0 SET lowestAutoCorrNum = MAX_AUTO_CORR SET K = 0 REPEAT WHILE K < number of BTS IF auto correlation number for BTS K > highestAutoCorrNum THEN SET highestAutoCorrNum = auto correlation number of BTS K SET highestIdx = K ELSE IF auto correlation number for BTS K < lowestAutoCorrNum THEN SET lowestAutoCorrNum = auto correlation number of BTS K SET lowestIdx = K END IF INCREMENT K END REPEAT DECREASE txGain for BTS highestIdx by TX_GAIN_STEP IF txGain for BTS highestIdx < MIN_TX_GAIN THEN SET txGain for BTS highestIdx = MIN_TX_GAIN END IF SET radio TX gain for BTS highestIdx = txGain for BTS highestIdx INCREASE txGain for BTS lowestIdx by TX_GAIN_STEP IF txGain for BTS lowestIdx > MAX_TX_GAIN THEN SET txGain for BTS lowestIdx = MAX_TX_GAIN END IF END

In summary, this method first determines the maximum and minimum auto-correlation values and records the indices for the corresponding BTS. Then, the transmit gain of the BTS with the highest auto correlation is reduced by TX_GAIN_STEP, and the transmit gain of the BTS with the lowest auto correlation is increased by TX_GAIN_STEP.

Finally, we show the performance of the transmit power balancing methods in practical outdoor propagation scenarios. The first scenario we considered is depicted in FIG. 42. Transmit power imbalance is caused by the two clients UE00, UE01 being in close proximity to BTS 10. The distribution of the condition number with and without power imbalance is shown in FIG. 43.

FIGS. 44-46 show the channel traces (SNDR, SSI and TSI) for three different cases: i) DIDO 2x2 without transmit power balancing method; ii) DIDO 2x2 with transmit power balancing method; iii) DIDO 4x2 with transmit power balancing method in combination with antenna selection method. The SSI threshold is set to −10 dB. We observe that transmit power balancing and antenna selection methods help improve the SSI trace and consequently the SNDR performance.

A different scenario is depicted in FIG. 47, where BTS 6 is the source of transmit power imbalance. FIG. 48 shows the condition number distribution with and without imbalance and FIGS. 49-51 depict channel traces for different algorithms as before. Once again, both transmit power balancing and antenna selection methods improve SSI and SNDR performance.

Another embodiment of the invention employs a combination of transmit power balancing and antenna selection algorithms. In this method, the extra antenna that provides the largest auto-correlation coefficient is removed and the conventional antenna selection algorithm is applied with the remaining extra antennas. For example, FIG. 52 shows the SER and CDF of the SSI for DIDO systems with 4 clients. The performance of DIDO 6x4 is significantly degraded when the system undergoes transmit power imbalance. By removing the BTS that yields a large auto-correlation coefficient, system performance is improved significantly as shown by the SER curve for DIDO 5x4.

Finally, FIG. 53 shows the statistics of the BTSs that have been selected for transmission and their relative usage. The last graph shows that DIDO 5x4 with transmit power balancing and antenna selection has zero usage of BTS1 because that BTS has been removed as a result of high auto-correlation value.

6. REFERENCES

-   [1] R. A. Monziano and T. W. Miller, Introduction to Adaptive     Arrays, New York: Wiley, 1980. -   [2] K. K. Wong, R. D. Murch, and K. B. Letaief, “A joint channel     diagonalization for multiuser MIMO antenna systems,” IEEE Trans.     Wireless Comm., vol. 2, pp. 773-786, July 2003; -   [3] M. Costa, “Writing on dirty paper,” IEEE Transactions on     Information Theory, Vol. 29, No. 3, Page(s): 439-441, May 1983. -   [4] U. Erez, S. Shamai (Shitz), and R. Zamir, “Capacity and     lattice-strategies for cancelling known interference,” Proceedings     of International Symposium on Information Theory, Honolulu, Hi.,     November 2000. -   [5] G. Caire and S. Shamai, “On the achievable throughput of a     multiantenna Gaussian broadcast channel,” IEEE Trans. Info. Th.,     vol. 49, pp. 1691-1706, July 2003. -   [6] M. Tomlinson, “New automatic equalizer employing modulo     arithmetic,” Electronics Letters, Page(s): 138-139, March 1971. -   [7] H. Miyakawa and H. Harashima, “A method of code conversion for     digital communication channels with intersymbol interference,”     Transactions of the Institute of Electronic -   [8] D.-S. Shiu, G. J. Foschini, M. J. Gans, and J. M. Kahn, “Fading     correlation and its effect on the capacity of multielement antenna     systems,” IEEE Trans. Comm., vol. 48, no. 3, pp. 502-513, March     2000. -   [9] J. P. Kermoal, L. Schumacher, K. I. Pedersen, P. E. Mogensen,     and F. Frederiksen, “A stochastic MIMO radio channel model with     experimental validation,” IEEE Jour. Select. Areas in Comm., vol.     20, no. 6, pp. 1211-1226, August 2002. -   [10] R. Chen, R. W. Heath, Jr., and J. G. Andrews, “Transmit     Selection Diversity for Unitary Precoded Multiuser Spatial     Multiplexing Systems with Linear Receivers,” IEEE Trans. on Signal     Processing, vol. 55, no. 3, pp. 1159-1171, March 2007. 

We claim:
 1. A method implemented within a multiuser (MU) multiple antenna system (MAS) comprising: communicatively coupling one or more centralized units to multiple distributed transceiver stations or antennas via a network, the network comprising wireline or wireless links or a combination of both; employing the network as a backhaul communication channel; and employing transmit antenna selection, user selection, or transmit power balancing to enhance channel space diversity within the MU-MAS; receiving wireless power concurrently by a plurality of user devices from the multiple antennas.
 2. The method as in claim 1 wherein a space selectivity indicator (SSI), time selectivity indicator (TSI), and/or frequency selectivity indicator (FSI) are used as measures of the channel diversity.
 3. The method as in claim 2 further comprising defining different subsets of transmit antennas and selects the subset that optimizes the SSI for transmission over the wireless links.
 4. The method as in claim 2 wherein a selection criterion is based on SSI, TSI and/or FSI thresholds derived from measurements.
 5. The method as in claim 2 wherein a minimum singular value of an effective channel matrix of every client device is used as SSI.
 6. The method as in claim 2 wherein a minimum singular value or the condition number of a composite channel matrix from all client devices is used as SSI.
 7. The method as in claim 2 wherein an absolute value of the sum of the complex channel gain from some or all transceiver stations is used as TSI.
 8. The method as in claim 2 wherein the SSI is used to measure and predict a shape in space around each of the plurality of user devices, each shape carrying independent and simultaneous non-interfering wireless power for the user device.
 9. The method as in claim 2 wherein an average SSI is used to select adaptively between fixed transmit antenna configuration and transmit antenna selection methods based on changing channel conditions.
 10. The method as in claim 2 wherein temporal correlation of the SSI is exploited to select an optimal antenna subset while reducing the computational complexity of the method.
 11. The method as in claim 3 wherein searching for the optimal antenna subset is suspended as soon as the first subset that satisfies the SSI threshold is found, thereby reducing computational complexity.
 12. The method as in claim 3 wherein only a limited number of antenna subsets are selected based on certain performance criterion as a means to reduce the computational complexity of the method.
 13. The method as in claim 1 wherein a base transceiver station (BTS) that overpowers the other BTSs reduces its transmit power to balance the power from all BTSs to the clients.
 14. The method as in claim 1 wherein a base transceiver station (BTS) that overpowers the other BTSs keeps its power level unaltered, the other BTSs increase their transmit power to balance the power from all BTSs to the client devices.
 15. The method as in claim 1 wherein a maximum auto-correlation coefficient of the covariance matrix is used as an indication of transmit power imbalance.
 16. The method as in claim 15 wherein thresholds of auto-correlation are defined to select between power balanced and imbalanced methods and those thresholds are obtained based on certain performance criterion.
 17. The method as in claim 16 wherein the selection is based on hysteresis loop and multiple thresholds of auto-correlation are defined for that hysteresis.
 18. The method as in claim 15 wherein the auto-correlation is mapped into a transmit gain value and that value is used to adjust the power of the BTSs.
 19. A multiuser (MU) multiple antenna system (MAS) comprising: one or more centralized units communicatively coupled to multiple distributed transceiver stations or antennas via a network; the network comprising wireline or wireless links or a combination of both, employed as a backhaul communication channel; and the MU-MAS employing transmit antenna selection, user selection and/or transmit power balancing are employed to enhance channel spatial diversity; receiving wireless power concurrently by a plurality of user devices from the multiple antennas.
 20. The system as in claim 19 wherein a space selectivity indicator (SSI), time selectivity indicator (TSI), and/or frequency selectivity indicator (FSI) are used as measures of the channel diversity.
 21. The system as in claim 20 further comprising defining different subsets of transmit antennas and selects the subset that optimizes the SSI for transmission over the wireless links.
 22. The system as in claim 20 wherein a selection criterion is based on SSI, TSI and/or FSI thresholds derived from measurements.
 23. The system as in claim 20 wherein a minimum singular value of an effective channel matrix of every client device is used as SSI.
 24. The system as in claim 20 wherein a minimum singular value or the condition number of a composite channel matrix from all client devices is used as SSI.
 25. The system as in claim 20 wherein an absolute value of the sum of the complex channel gain from some or all transceiver stations is used as TSI.
 26. The system as in claim 20 wherein the SSI is used to measure and predict a shape in space around each of the plurality of user devices, each shape carrying independent and simultaneous non-interfering wireless power for the user device.
 27. The system as in claim 20 wherein an average SSI is used to select adaptively between fixed transmit antenna configuration and transmit antenna selection systems based on changing channel conditions.
 28. The system as in claim 20 wherein temporal correlation of the SSI is exploited to select an optimal antenna subset while reducing the computational complexity of the system.
 29. The system as in claim 21 wherein searching for the optimal antenna subset is suspended as soon as the first subset that satisfies the SSI threshold is found, thereby reducing computational complexity.
 30. The system as in claim 21 wherein only a limited number of antenna subsets are selected based on certain performance criterion as a means to reduce the computational complexity of the system.
 31. The system as in claim 19 wherein a base transceiver station (BTS) that overpowers the other BTSs reduces its transmit power to balance the power from all BTSs to the client devices.
 32. The system as in claim 19 wherein a base transceiver station (BTS) that overpowers the other BTSs keeps its power level unaltered, the other BTSs increase their transmit power to balance the power from all BTSs to the clients.
 33. The system as in claim 19 wherein a maximum auto-correlation coefficient of the covariance matrix is used as an indication of transmit power imbalance.
 34. The system as in claim 33 wherein thresholds of auto-correlation are defined to select between power balanced and imbalanced systems and those thresholds are obtained based on certain performance criterion.
 35. The system as in claim 34 wherein the selection is based on hysteresis loop and multiple thresholds of auto-correlation are defined for that hysteresis.
 36. The system as in claim 35 wherein the auto-correlation is mapped into a transmit gain value and that value is used to adjust the power of the BTSs. 